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We investigate the instabilities and bifurcations of traveling pulses in a model excitable medium; 
in particular we discuss three different scenarios for the loss of stability resp. the disappearance 
of stable pulses. In numerical simulations beyond the instabilities we observe replication of pulses 
("backfiring") resulting in complex periodic or spatiotemporally chaotic dynamics as well as modu- 
lated traveling pulses. We approximate the linear stability of traveling pulses through computations 
in a finite albeit large domain with periodic boundary conditions. The critical eigenmodes at the 
onset of the instabilities are related to the resulting spatiotemporal dynamics and "act" upon the 
back of the pulses. The first scenario has been analyzed earlier [M. G. Zimmermann et al, Physica 
D 110, 92 (1997)] for high excitability resp. low excitation threshold: it involves the collision of a 
stable pulse branch with an unstable pulse branch in a so called T-point. In the frame of traveling 
wave ordinary differential equations, pulses correspond to homoclinic orbits and the T-point to a 
double heteroclinic loop. We investigate this transition for a pulse in a domain with finite length and 
periodic boundary conditions. Numerical evidence of the proximity of the infinite-domain T-point 
in this setup appears in the form of two saddle-node bifurcations. Alternatively, for intermediate 
excitation threshold, an entire cascade of saddle-nodes causing a "spiraling" of the pulse branch 
appears near the parameter values corresponding to the infinite domain T-point. Backfiring ap- 
pears at the first saddle-node bifurcation, which limits the existence region of stable pulses. The 
third case found in the model for large excitation threshold is an oscillatory instability giving rise 
to "breathing", traveling pulses which periodically vary in width and speed. 

PACS numbers: 03.40.Kf, 47.54. -fr, 82.40. Ck 



I. INTRODUCTION 

One-dimensional excitable media exhibit nonlinear 
traveling waves such as wave trains and solitary pulses. 
Examples include concentration waves in chemical reac- 
tions in solution and on surfaces , signal propagation 
in neurons and in cardiac tissue ^ . These processes are 
often modelled by reaction-diffusion (RD) equations of 
activator-inhibitor type |^,|^ . A pulse is a localized struc- 
ture; it may result from a finite amplitude perturbation 
of a linearly stable rest state. The pulse shape usually de- 
cays exponentially as a function of the distance from the 
pulse center. Pulses can be analytically approximated in 
the asymptotic limit where the dynamics of the activator 
is much faster than that of the inhibitor variable, and 
where the activator variable diffuses, but the inhibitor 
does not. There, it can be shown that a pulse exists and 
IS stable [|. 

More recently, spatiotemporal chaos has been found 
in a variety of excitable model systems. Several exam- 
ples have been reported in models whose kinetics possess 
three distinct homogeneous steady states (fixed points); 
apart from the stable rest state, these media also exhibit 
two unstable fixed points |^-^. One example arises in a 
model for catalytic CO oxidation, where the inhibitor is 



a so-called surface structure variable; its kinetic nuUcline 
displays a sigmoidal shape - threshold dependence on the 
activator variable (adsorbate coverage) In contrast to 
the standard linear dependence of the inhibitor nuUcline 
on the activator, this functional form leads to the addi- 
tional fixed points mentioned above. Similar behavior is 
often seen in models in physiology describing the dynam- 
ics of the membrane voltage (activator) controlled by so 
called gating variables for the ion channels (inhibitor). 
The kinetics of these gating variables display again a 
threshold- type, sigmoidal dependence on the membrane 
voltage, often leading to three (spatially homogeneous) 
fixed points for the PDE Q . If only a single gating vari- 
able is involved, these models are qualitatively similar 
to the equations for catalytic CO oxidation studied here; 
a good example is provided by the Morris-Lecar model 
used to describe the membrane potential in a barnacle 
muscle fiber fl0[| . 

If the inhibitor kinetics are fast enough, the CO oxi- 
dation model displays an instability that has been col- 
loquially named backfiring of pulses Related phe- 
nomena include the wave-induced chemical chaos found 
in the Gray-Scott model for an autocatalytic chemical 
reaction ||^, as well as complex behavior in amplitude 
equations describing dynamics near a Takens-Bogdanov 
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(TB) point In addition, even reaction-diffusion media 
whose kinetics are cliaracterized by a single fixed point, 
sometimes display pulse instabilities and backfiring under 
excitable conditions |ll|Jf 2t] . An analysis of the travel- 
ing wave ordinary differential equations derived from the 
original PDEs for the CO oxidation model reveals the ba- 
sic mechanism for the destruction of stable pulses leading 
to complex behavior: it involves a so-called T-point |13| ] 
as well as a rich web of bifurcations of pulse solutions 
10. More recently, a similar bifurcation structure has 
been found in a model for intracellular calcium waves in 
pancreatic acinar cells p5| and in the ultrarcfractory ver- 
sion of the FitzHugh-Nagumo model The resulting 
complex dynamics is often governed by a coherent struc- 
ture described as a "wave emitting front" [l^ , a nonlinear 
front involving a spatially uniform unstable state that in- 
vades a spatially uniform stable one. The unstable state 
behind the propagating front evolves into spatiotemporal 
chaos. Similar behavior has been seen in the amplitude 
equations in the neighborhood of the TB point, and was 
named "chaotic nucleation" by those authors . 

In this paper we investigate instabilities and bifurca- 
tions in a model of a catalytic surface reaction Us- 
ing a computer-assisted approach we calculate solution 
branches in large, finite domains, under periodic bound- 
ary conditions (which we will, with slight abuse of ter- 
minology, refer to as pulses) and perform linear stability 
analysis of these pulse solutions. The eigenvalues of the 
corresponding Jacobian reflect the growth exponents of 
perturbation modes. The emphasis of the present paper 
is on the form of the spectra and in particular on the 
shape of the critical modes; this is in contrast to ear- 
lier work which focused mostly on the pulse bifurcation 
structure (l3[|l|]. 

A homogeneous steady state of reaction-diffusion 
equations in an infinite system possesses a continuous 
spectrum, reflecting the growth rates of perturbation 
modes. The eigenmodes are harmonic functions, whose 
wavenumber parametrizes the corresponding eigenvalues. 
In a very long or infinite system, a single pulse can con- 
ceivably be considered (in the appropriate norm) as a 
perturbation of this uniform solution. Under these cir- 
cumstances, the modes far away from the pulse are still 
delocalized, harmonic wave modes; the continuous spec- 
trum remains essentially unchanged. Furthermore, ad- 
ditional perturbation modes exist that are localized at 
the site of the pulse and decay exponentially away from 
it. The corresponding eigenvalues are discrete, in con- 
trast to the continuous band of eigenvalues belonging to 
the nonlocalized modes. In infinite or periodic systems, 
the Goldstone mode is an example of such a localized 
mode. In a spatially homogeneous system it is given by 
the spatial derivative of the pulse profile, and accounts 
for a shift in space (due to translational invariance); the 
corresponding eigenvalue is zero. 

A bifurcation or instability of a given solution upon 
change of a single control parameter is typically accom- 
panied by a single real eigenvalue or a pair of complex 



conjugate eigenvalues crossing the imaginary axis. For 
a real eigenvalue, a saddle-node bifurcation is the most 
common case, while alternative bifurcations like the tran- 
scritical or pitchfork bifurcations require certain symme- 
tries for the equations and the solution. The case of a 
drift pitchfork bifurcation of a stationary pulse was stud- 
ied in p^ . A destabilizing perturbation grows in an os- 
cillatory manner if the eigenvalues are complex and in a 
non-oscillatory manner if the corresponding eigenvalue is 
real. The first case leads to a Hopf bifurcation, while the 
latter case may correspond either to a saddlc-node or to 
a transcritical or pitchfork bifurcation. 

Three examples of dynamics where localized initial 
conditions beyond the uniform solution excitation thresh- 
old do not eventually evolve into stable pulses, are shown 
in Fig. ^ For high excitability resp. small excitation 
threshold, pulse-like initial conditions evolve in a travel- 
ing pulse that splits off new pulse-like excitations trav- 
eling in the opposite direction. This is the phenomenon 
termed backfiring in It can eventually evolve in spa- 
tiotemporally chaotic (I) or complex periodic (II) fash- 
ion. For large excitation threshold, modulated traveling 
pulses are also seen (III) . Their shape undergoes a breath- 
ing, periodic variation. 




FIG. 1. Space-time plots from numerical integration of 
Eqs. (1) showing the time evolution of pulses at parame- 
ter values beyond the instability onset. A stable pulse solu- 
tion for a subcritical parameter value, characterized by rest 
state A, was used as the initial condition. I: Backfiring in 
the immediate vicinity of a T-point. The resulting behavior 
is non-periodic for the given system length and initial con- 
ditions. 11: Backfiring after a saddle-node bifurcation. The 
resulting behavior is periodic in time for the given system 
length and initial conditions. Ill: A supercritical Hopf bifur- 
cation leads to modulated traveling waves. The pulse shows 
a periodic oscillation of shape and speed; variations appear 
mainly at its trailing edge. Black denotes high values of 
It, white corresponds to it = 0. Parameters: I: fe = 0.07, 
e = 0.1075,i = 100, AT = 119.2, II: b = 0.15, e = 0.0931, 
L = 100, AT = 238.42, III: b = 0.2, e = 0.062, , L = 50, total 
integration time AT = 238.42. 
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II. MODEL AND METHODS 

We investigate a model of activator-inhibitor type orig- 
inally derived for CO oxidation on Pt(llO) in a pa- 
rameter regime where the kinetics give rise to one stable 
and two unstable steady homogeneous solutions. This 
model is related to the FitzHugh-Nagumo system and 
describes the interaction of a fast activator u and a slow 
inhibitor variable v: 

Otu — of.u + -u[l - u)(u j , 

e a 




(1) 



< u < 1/3 
1)2, 1/3 < u < 1, 

1 < u 



with X € [0, L] and periodic boundary conditions. 

The time scales ratio e > is used as the control pa- 
rameter. The case e — > corresponds to the aforemen- 
tioned asymptotic limit where stable pulses are expected. 
The parameter b controls the excitability threshold of the 
system: the value of b is proportional to the magnitude of 
the critical perturbation which will trigger a pulse. The 
value of a is fixed at 0.84 throughout the paper. We will 
vary b, thus varying the excitability. In the parameter 
range considered here, three relevant fixed points exist: 
the stable state A — (0,0), the saddle B = (6/a,0), and 
unstable focus C. 

We shall now examine the stability of pulses travel- 
ing in a background of the stable rest state A. Since 
we investigate solutions moving with fixed velocity c, the 
analysis of their stability is performed in the comoving 
frame z :— x — ct: 



1 



dtu — d^u + cdzU H — u{l — u){u - 
dtv = cdzV + f{u) 



e 

V . 



(2) 



In this frame, traveling waves with speed c correspond 
to time independent, steady solutions. Linearization of 
these equations around a stationary solution uq{z), Vo{z) 
yields an eigenvalue problem for small perturbations 

(r(z,i),s(z,i))«(r-(z),s(z))e^*: 
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d'i + cd, + g,{z) 
duf (wo) 
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92{z) 
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with 



gi{z) = -- \ uq{uq - 1) + (uo - - 

e \ a 

92[z) = — (uo - 1). 
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The linear stability problem amounts to the determina- 
tion of the spectrum of the Jacobian M in (0). For the 



homogeneous steady states A or B, at least one of the 
off-diagonal matrix elements is zero. Thus, the diagonal 
elements of the matrix M suffice to compute the spec- 
trum. For the A steady state it is: 



\ai ^ + ick, 2 = — 1 + ick, (4) 

ea 

and for the B steady state it is: 

Xb 1 ~ I 1 I — fc^ + ick, As 2 = — 1 + ick, 

ea \ a J 



(5) 



where k is the wavenumber of the perturbation. In 
the case of periodic boundary conditions studied here, 
k — applies. The real part of the eigenvalues A2 is 



— 1. The eigenvectors are then {r,s)'^ — (1,0) 
(0,l)^e*^■^ 

In general, traveling wave solutions uq{z),vo{z) and 
the eigenfunctions of the Jacobian Ai are not available 
in closed form. Thus, the problem has to be approxi- 
mated numerically. We approximate pulses by comput- 
ing steady solutions in a finite length system and the trav- 
eling frame through a pseudospectral discretization of 
Eqs. ^ with periodic boundary conditions and Newton- 
Raphson iterations. The velocity c, which is not known 
a priori, is formally an additional variable along with the 
Fourier coefficients of the solution. One additional pin- 
ning condition singles out one of the infinitely many solu- 
tions existing due to translational invariance and allows 
the numerical computation of the speed. The eigenfunc- 
tions and the spectrum of the Jacobian are obtained in 
Fourier space resulting from a 200 {resp. 250 for case 
I) mode decomposition of the stationary solution. The 
zero eigenvalue, which always exists due to the transla- 
tional symmetry of the problem, is used as a numerical 
accuracy check and has been obtained with a precision of 
lO^'^ or better. The time evolution of unstable solutions 
was computed using an explicit finite difference scheme 
to solve Eqs. (|l|), discretizing space in 1024 points and 
using a time step of At — 0.0122. 

Because traveling solutions fulfill the conditions dtU = 
dfV = in the comoving frame, they can also be ob- 
tained in the traveling wave ordinary differential equa- 
tions (TWODE) following from Eqs. (|): 



du 

dz ' 

dw 1 

— = -CW + -U{U - 1)(U - Uthj 



(6) 



dv 
dz 



= {v~f{u))/c, 



with uth = {b + v)/a. In this frame, a homogeneous solu- 
tion corresponds to a fixed point, a pulse to a homoclinic 
orbit, and a front to a heteroclinic orbit. Consequently, 
in the parameter range studied here, three relevant fixed 
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points exist: A 
C. 



(0,0,0), B = (Va,0,0), and the focus 
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case III {b = 0.20). We proceed as follows: first, branches 
of pulse solutions on a ring, representative of ,,true" pulse 
solutions in a infinite domain, are presented for the rel- 
evant range of the control parameter e; we characterize 
these pulses by their speed c (Fig. |2[ ||b). We then show 
pulse profiles and spectra at selected values of e shortly 
before and after the onset of instability (Figs. MiHi 
and present the destabilizing mode (r, s)'^ (Figs.^, 0, i ). 
Representative post-instability spatiotemporal dynamics 
can be found in Fig. |l]. 

Fig. H shows pulse speed as a function of the parameter 
e for the three cases. The thick lines correspond to stable 
pulses with background state A, while dashed lines corre- 
spond to unstable pulses with background state B. The 
transitionpoint between the two families is the so called 
T-point [P3| , denoted by a T. This point is defined as a 
double hcteroclinic connection between the fixed points 
A and B in the frame of the traveling wave ODEs, Eqs. 
(^. In the frame of the original equation (^, these hcte- 
roclinic orbits correspond to fronts. The branch of pulses 
corresponding to homoclinic orbits to A in the TWODE 
may (Fig. |2[ll. III) or may not (Fig. ^.I) spiral into the 
T-point. Spiraling is observed when two of the eigenval- 
ues of the linearization of the TWODE around the fixed 
point B are complex conjugate [p^ . 
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FIG. 2. Pulse speed as a function of the parameter e. I: 
b = 0.07, II: b = 0.15, and III: b = 0.2. Thick lines are pulses 
with rest state A (i.e. pulses in a large box with periodic 
boundaries representative of pulses to the state A on the in- 
finite line corresponding to the homoclinic orbits to A); thin 
lines are pulses with rest state B. Solid (dashed) lines denote 
stable (unstable) branches. A T marks our approximation of 
the double heteroclinic connection point (the T-point) where 
pulses with rest state A "collide" with pulses with rest state 
B. H denotes Hopf bifurcations, SN denotes saddle-node bi- 
furcations. In the cases II and III, the branch of pulses with 
rest state A spirals into the T-point. This is not the case for 
I, compare Fig. Hb. 



III. RESULTS 

We consider three cases, at increasing values of the ex- 
citation threshold, controlled through the parameter b at 
fixed a = 0.84: Case I (6 = 0.07), case II {b= 0.15) and 



A. Case I 



Fig. H.I shows speed as a function of the parameter e 
for b = 0.07 computed for pulses in a large system with 
periodic boundaries. The picture appears at first glance 
identical with the result found in a continuation of homo- 
clinic orbits in the TWODE The branch of pulses 
to A (thick line) does not spiral into the T-point. To 
be exact, the T-point only exists for pulses on the infi- 
nite line; what - at the resolution of our picture - still 
appears as a T-point will be discussed in more detail be- 
low. To understand the nature of the instability in this 
case, we focus on the solutions near the T-point value 
CT « 0.10744. 

As mentioned above, for an infinite system, the T-point 
corresponds to a double heteroclinic connection in the 
traveling wave ODEs, Eqs. (||). Close to ex the orbits 
homoclinic to A also approach the fixed point B (and 
vice versa) ; this approach to B makes the dynamics in z 
progressively slower. Thus, close to et, pulses with rest 
state A (resp. B) will also exhibit extensive regions close 
to B (resp. A) , practically giving rise to fronts between A 
and B within the pulse profile. When we perform a con- 
tinuation for a pulse in a ring of finite length, the space 
on the ring is apart from the excitation plateau divided 
between residence close to A and residence close to B, so 
that the total period is constant as we vary e. 
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FIG. 3. a) Four pulse solutions on the ring from the mid- 
dle branch in (b) illustrating the gradual transformation from 
solutions identifiable with pulses with rest state A {u — 0) to 
unstable solutions identifiable with pulses with rest state B 
{u = 0.083). The B domain in the waveform becomes wider 
upon decrease of e, along the middle branch in (b). b) Bifurca- 
tion diagram, with respect to e, of pulses for case I, exhibiting 
two saddle-nodes. The thick solid line can be associated with 
stable pulses homoclinic to A; the thin dashed line can be 
associated with unstable pulses homoclinic to B. The thick 
dashed line corresponds to the transition region between these 
two cases; it constitutes the incarnation (for our finite, large 
ringlength continuation) of the T-point. When the upper sad- 
dle-node bifurcation (SNi) takes place, our finite ringlength 
solution could be still described as an approximation of a pulse 
to A but with a small B shoulder; the converse description 
holds at the lower saddle-node bifurcation {SN2). It is reason- 
able to consider as most representative of the infinite-domain 
T-point the location, along this thick-dashed line, where the 
pulse solution contains comparable large patches close to A 
and close to _B - roughly the middle of this middle-branch. 

To study this behavior, numerical continuation tech- 
niques are needed. We computed the stationary solu- 
tions of Eqs. on a ring of length L = 200 with 250 
modes, 1024 collocation points and a parameter step size 
Ae — 0(1G~^°). Periodic approximations of pulse so- 
lutions with rest state A are shown in Fig. |^a as they 
approach the T-point and undergo the gradual transi- 
tion to pulse solutions with rest state B. One can clearly 
recognize the increasing domain B — (b/a, 0) at the back 
of the pulse. As the state B is unstable, one might expect 
that the pulse in the ring loses stability as soon as the B 
domain gets large enough. However, one should be aware 
of the fact that the B domain is moving with the speed of 
the pulse. It is possible that the B domain is only convec- 
tively unstable in the comoving frame. In other words. 



the perturbations growing on the B state in a stationary 
frame may spread slower than the pulse motion, and the 
pulse in an infinite domain will be stable. A mathemat- 
ically precise description of this phenomenon has been 
given by Sandstede and Scheel . A related result has 
been obtained by Nii |2^] , who shows that eigenvalues ac- 
cumulate in the area bounded by the essential spectra of 
A and B. The opposite case, i.e. perturbations on the B 
plateau spread faster than the pulse speed, is described 
in the next subsection (case II). 

Typical spectra are shown in Fig. ^ for solutions along 
the c — e-branch. To facilitate a better comparison, the 
continuous spectra of the rest states A (Eqs. (||)) and B 
(Eqs. (^) are depicted as solid lines; the computed eigen- 
values are denoted by circles. Fig. ^ shows the spectrum 
of a solution approximating a stable pulse with rest state 
A, while Fig. Qf shows that of a solution approximating 
an unstable pulse with rest state B. In both cases, the 
parameter e is "far" from the T-point bifurcation and the 
eigenvalues belonging to nonlocalized eigenvectors com- 
pare well with the continuous spectrum of the respective 
rest state. Figs, ^-e show the gradual transition between 
the two cases. Fig. ^ shows the spectrum for a stable 
solution for which the domain B is just starting to appear 
in the back of the pulse. Its width is Lb ^ 5. Several 
discrete eigenvalues have appeared on the negative real 
axis; a similar spectrum has been predicted and found by 
Sandstede and Scheel |l9| ]. 

The details of the spectrum tranformation are as fol- 
lows: Initially (far to the left of the T-point in Fig. |^), 
the pulse spectrum contains two discrete complex conju- 
gate eigenvalues (which remain more or less unchanged 
through the bifurcation diagram) as well as two "pieces" 
of continuous spectrum: a parabola, and a vertical line. 
For our calculations in the finite domain, initially the 
parabola does have a "tip" (an eigenvalue on the real 
axis) while the vertical line does not have an eigenvalue on 
the real axis. As we vary e the parabola approaches the 
straight line; eventually the parabola appears as cross- 
ing the line. While this occurs, a number of complex 
conjugate pairs of eigenvalues emerge from the parabola 
and line approximating the essential spectrum. Fig. ^ 
shows the situation where the first of these pairs has ap- 
peared. Upon further increase of e, these pairs collapse 
on the negative real axis where they split and become 
real. These real eigenvalues lie in the absolute spectrum 
of state B (see 1 19 ) ; they appear when the solution has a 
visible B "shoulder" behind the pulse, see Fig. Ha. As e 
varies these real eigenvalues become complex conjugate. 
Several complex eigenvalue pairs of the parabola give the 
appearance of a secondary parabola - this is the first inti- 
mation of what will become the essential spectrum of B. 
Only two of the "momentarily real" eigenvalues remain 
real - one of them forms the tip of the "new parabola" , 
while the other forms the tip of the "old parabola" . 

Gradually, as the B plateau grows, two distinct phases 
of spectrum movement are observed. First this new B- 
associated parabola moves to the right in the complex 
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plane and at some point it starts crossing the imaginary 
axis. That is precisely the first saddle-node bifurcation 
we observe - the critical eigenvalue is the tip of this sec- 
ondary parabola, whose origin we just discussed, see Fig. 

Continuing further on the middle, already unstable 
branch, the _B-domain gets wider and many secondary 
(apparently Hopf) bifurcations occur as the eigenvalues 
of this "discretized" parabola cross successively the imag- 
inary axis. When both A- and i?-plateaus are equally 
present, one expects to see echoes of both A- and B- 
continuous spectra in the solution spectrum, and that is 
indeed seen in Fig. Near the lower saddle-node bi- 
furcation, the i3-plateau is almost as wide as the system 
length. There, one may describe the solution more rea- 
sonably as an approximation to a pulse with rest state 
B with a short A domain at its front. The eigenvalue 
spectrum at this situation is depicted in Fig. The sec- 
ond (lower) turning point occurs when the "discretized" 
parabola of eigenvalues that corresponded initially to 
continuous spectrum of A moves to the right and the real 
eigenvalue at the tip of this discretized parabola crosses 
the imaginary axis. 
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the final pure-B parabola has no real 
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FIG. 4. Spectra of pulses in case I for different values of 
e. The solid lines correspond to the continuous spectrum of 
the rest state A, Eqs. H (left parabola) resp. rest state B 
(right parabola), Eqs. ^ The pictures show the typical change 
of the spectrum along the solution branch c-e. Parameters: 
a = 0.84, b = 0.07, L = 200 (L = 400 for picture d), a) 
e = 0.106425824, b) e = 0.107446965, c) e = 0.107446995547, 
d) e = 0.1079, e) e = 0.107446995527, f) e = 0.109936154. 



Eventually, the old A-parabola eigenvalues merge with 
the i?-parabola as the A-shoulder gradually disappears. 
Fig. |f shows the spectrum for a solution approximating 
a pulse with rest state B far from the T-point bifurcation 
for a high values of e. In the end, the two "tips" (the old 
and the new one) collide and become again a complex 



Sandstede and Scheel in |T^] have proven that the spec- 
trum of pulses with rest state A with a "shoulder" of B 
state as in Fig. ||a is comprised of the essential spectrum 
of the homogeneous background state A plus eigenval- 
ues lying in the absolute spectrum Q of B located on 
the real axis and bounded by a maximum value which 
is negative. We deviate from their picture, when the B 
shoulder becomes relatively large. In other words, their 
result on the infinite line is reproduced by our numerical 
stability analysis in the ring, if the condition Lb « L 
holds. The spectrum transformations, however, reflect 
the influence of the periodic boundary conditions and 
the finite length. Because of the violation of the above 
assumption near the T-point (see Fig. 3a) the spectra we 
find for a large box are not well described as the union 
of the essential spectrum of A and the absolute spectrum 
of B. The results of jl^ also imply that if one considers 
pulses on the infinite line instead of the large wavelength 
approximation with periodic boundary conditions used 
here, one will find a monotonic approach of the pulse 
branches to the T-point. In their scenario no eigenvalues 
are crossing the real axis upon approaching the T-point 
(in contrast to the computations on the ring) and no bi- 
furcation (like the saddle-node bifurcation found for the 
ring) can occur. 



Focusing on the first instability {SNi), the critical 
eigenmode corresponding to the saddle-node appears lo- 
calized at the back of the pulse (see Fig. 5). In addition, 
it is acting mainly on the activator (Fig. 5b); we ob- 
serve that new pulses can now split off the existing pulse 
("backfiring"). Fig. || shows the time evolution beyond 
the saddle-node bifurcation in a comoving frame. The 
initial condition is the stable pulse just before the bifur- 
cation. Fig. I^a shows how a perturbation with support 
over the i3-plateau grows with time, while the shape of 
the preceding front essentially does not change. The back 
(what starts as the i?-plateau) grows until its maximum 
reaches the excited state u = 1. The excited domain 
widens, but its plateau is unstable against oscillating per- 
turbations Therefore, the situation rapidly evolves, 
and the growing excited domain breaks down, giving rise 
to two new seeds for ^-pulse like entities. This behav- 
ior has been described by the term "backfiring" []6|,p^. 
Backfiring occurs repeatedly, and the newly generated 
pulse-like structures annihilate upon collision with simi- 
lar objects traveling in the opposite direction; this inter- 
play of instability, new pulse generation and annihilation 
gives rise to the non-periodic space-time behavior shown 
in Fig. II. 
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FIG. 5. a) Unstable pulse solution in case I in the 
vicinity of the saddle-node bifurcation SNi, Fig. ^b. b) 
Destabilizing eigenmode. Parameters: a — 0.84, b — 0.07, 
e = 0.1074469956, L = 200. Only part of the domain is 
shown. 
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FIG. 6. Time evolution for a value of e after the sad- 
dle-node bifurcation SNi in a comoving frame. The initial 
condition is the stable pulse-like solution just before the bi- 
furcation (thick black line in a)). The thin lines in a) demon- 
strate how the shape of the solution changes only at the back 
of the pulse-like solution. A localized perturbation grows in 
amplitude and width in the course of time and approaches 
the rest state C, where w = 1. Further evolution of the dy- 
namics is given by the dashed line in picture b). Due to 
oscillatory instability of the homogeneous state C, the corre- 
sponding domain cannot become large. A breakdown leads 
to the generation of two further pulse-like states traveling in 
opposite directions (thin solid line). 

Altogether, our results give an interesting and well 
resolved picture of stability of pulses on a large finite 
ring near the T-point. In contrast, earlier computations 
|13| , P^ have used coarser steps in the continuation algo- 
rithms and smaller domains. The present computations 
show how the pulse solutions change in a gradual fash- 
ion; and that two distinct solutions branches show an 
extremely narrow hysteresis; the middle branch medi- 
ates the transition between the two different pulse types. 
The main changes appear in an extremely tiny region 
of parameter space and can therefore only be captured 
with careful numerics. The spectrum transformations 
near parameters that exhibit a T-point in the correspond- 
ing TWODE appear in a similar fashion for case II and 
III studied below. There, however, the T-point is not 
involved in the disappearance of stable traveling pulses 
through instability and bifurcation. Thus, we will not 
discuss it in these cases. 



B. Case II 

Fig. |[ll shows the c — e diagram for b = 0.15. As 
in case I, we again observe here the transformation from 
pulses with rest state A to pulses with rest state B at 
what appears like a T-point. The branch corresponding 
to "pulses to A" spirals into this T-point. This behavior 
is caused by imaginary eigenvalues of the fixed point B 
in the TWODE at the T-point conditions, and has been 
predicted from general arguments . Thus, the branch 
of initially stable pulses with rest state A undergoes a 
sequence of saddle-node bifurcations upon approaching 
the T-point. This sequence of saddle-node bifurcation is 
consistent with the results of Sandstede and Scheel . 
The spectrum of A-pulses near the T-point in the infinite 
system approaches the union of the essential spectrum of 
A and the absolute spectrum of B. In this case, how- 
ever, the absolute spectrum of B contains real positive 
eigenvalues. One may, therefore expect infinitely many 
saddle-node bifurcations as the T-point is approached. 
In the stability-relevant first saddle-node bifurcation, the 
solution branch of stable wavetrains with rest state A 
turns around and becomes unstable. As is required for a 
saddle-node bifurcation, a single real eigenvalue crosses 
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the imaginary axis. This can be seen in Fig. left col- 
umn. The two spectra correspond to pulse-like solutions 
before and after the saddle-node bifurcation along the so- 
lution branch. The unstable pulse and the destabilizing 
eigenmode are shown in Fig. ^, right column. Once more, 
the destabilizing mode affects primarily the back of the 
pulse. It is worth noting, that upon further continuation 
the branch of pulses to A spirals towards the T-point 
related situation over a cascade of saddle-node bifurca- 
tions. Each saddle-node bifurcations adds an additional 
positive eigenvalue to the spectrum. This is in contrast 
to the study in |Q , where a single eigenvalue repeatedly 
crosses zero along a spiraling pulse branch with a cascade 
of saddle-node bifurcation. 
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FIG. 7. Case II: a) Eigenvalues at e = 0.0927 before a 
saddle-node bifurcation of a pulse decaying into the rest state 
A; b) eigenvalues at e = 0.09299 after this bifurcation. A dis- 
crete eigenvalue crosses the imaginary axis, c) Unstable solu- 
tion and d) corresponding destabilizing eigenmode (r, 5*)"^ at 
e = 0.09299. Parameter: a = 0.84, b = 0.15, L = 100. 

Numerical simulations of the model Eqs. ^ shortly af- 
ter the saddle-node bifurcation, for values of e for which 
no solution with rest state A exists, exhibits the phe- 
nomenon we termed above "backfiring" . In the tran- 
sients, our pulse-like object generates near its back other 
pulse-like entities traveling in the opposite direction, see 
Fig. |l|.II. After this transient period, the resulting spa- 
tiotemporal pattern in our finite domain becomes peri- 
odic in time. Simulations show, though, that this ob- 
servation depends on the initial conditions: non-periodic 
patterns like the one shown in Fig. |l|.I may also appear 
for the same parameter values. 



C. Case III 



Fig. 0.III shows the c — e for & = 0.2. Here, the pulse 
solution with rest state A becomes unstable through a 
Hopf bifurcation even before the first saddle-node bifur- 
cation is reached. The eigenvalue spectrum is shown in 
Fig. H, left column, on both sides of this Hopf bifurcation. 
It can be seen that one discrete pair of complex conju- 
gate eigenvalues crosses the imaginary axis. The second 
column of Fig. || shows the pulse-like solution after the 



bifurcation as well as the real and imaginary parts of the 
critical eigenmode (r, s)^. Note that the perturbation 
due to the eigenmode is localized and has its main con- 
tribution - once more - at the back of the pulse. As the 
corresponding eigenvalues are complex conjugate, the re- 
sulting perturbation oscillates with time at a frequency 
given by the imaginary part of the eigenvalues. 
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FIG. 8. Case III: a) Eigenvalues at e = 0.053 before a Hopf 
bifurcation; b) eigenvalues at e = 0.072 after this bifurcation. 
A discrete pair of complex conjugated eigenvalues crosses the 
imaginary axis, c) Unstable pulse-like solution as well as d) 
real and e) imaginary part of the corresponding destabiliz- 
ing eigenmode (r, s)"^ at e = 0.065. Parameters: a — 0.84, 
b = 0.2, L = 50 

Numerical integration of the model Eqs. (0) shows that 
the Hopf bifurcation which leads to destabilization is su- 
percritical iQ. The resulting pattern after destabiliza- 
tion consists of a modulated traveling pulse which oscil- 
lates in time, especially at its back (compare the simula- 
tion shown in Fig. ^III.) 



IV. CONCLUSION 

We have investigated the transition from stable pulse 
propagation to various regimes of more complicated spa- 
tiotemporal dynamics, namely modulated pulses, peri- 
odic and chaotic pulse backfiring. In all three cases, the 
transition can be explained by either a Hopf instability 
(modulated pulses) or a saddle-node bifurcation (leading 
to backfiring) of the stable pulse solution. In all cases, 
the transition is connected with either a single or a pair 
of complex conjugate discrete eigenvalue(s) with zero real 
part(s). In a finite domain with periodic boundary con- 
ditions - the typical experimental setup for investigation 
of pulses - spectra change continuously near the T-point 
in the fashion described in case I. The form of the corre- 
sponding critical eigenmode(s) allows some insight into 
how pulses become unstable resp. evolve in space and 
time. The dynamics in general still contains mostly prop- 
agating localized pulse-like structures whose evolution is 
governed by the unstable eigenmode (s) in the Hopf case 
(modulated) or by the critical eigenmode of the saddle- 
node bifurcation. Typically, the critical eigenmodes have 
support at the back of the pulse. 
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The results here should carry over to models with sim- 
ilar phenomenology mentioned in the introduction. Pre- 
liminary results p5| show that the transition to wave- 
induced chemical chaos in the Gray-Scott model Q is 
also caused by a saddle-node bifurcation of pulses near a 
T-point. T-points can only be found in systems with mul- 
tiple homogeneous fixed points {e.g. one stable rest state 
and two additional unstable steady states). The com- 
plex behavior seen in media with a single stable fixed 
point ||ll|,|l2| may be caused through a different mech- 
anism. For a model of the catalytic NO-CO reaction, 
upon change of control parameter first modulated travel- 
ing waves are seen and then periodic backfiring is found. 
This cannot be due to a T-point, but may instead be 
caused by a global bifurcation of the periodic modulated 
pulses - a scenario already suggested in a study of the 
present model with different control parameters . Fi- 
nally, it is important to note that a simple instability of 
a finite wavelength pulse solution, like the Hopf bifurca- 
tion in case III, only leads to a modulation of the shape, 
while a saddle-node bifurcation limits the existence of a 
certain type of pulses and may give rise to replication 
of pulse- like structures (cases I, II). The role of saddle- 
node bifurcations in the replication of pulses and creation 
of space-time defects has recently been investigated in 
the Gray-Scott model and in the Complex-Ginzburg- 
Landau equation [ p7| . 
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